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In this paper, an exact unitary transformation is examined that allows for the construction 
of solutions of coupled nonlinear Schrodinger equations with additional linear field coupling, from 
solutions of the problem where this linear coupling is absent. The most general case where the 
transformation is applicable is identified. We then focus on the most important special case, namely 
the well-known Manakov system, which is known to be relevant for applications in Bose-Einstein 
condensates consisting of different hyperfine states of 87 Rb. In essence, the transformation con- 
stitutes a distributed, nonlinear as well as multi-component generalization of the Rabi oscillations 
between two- level atomic systems. It is used here to derive a host of periodic and quasi-periodic 
solutions including temporally oscillating domain walls and spiral waves. 



Introduction. The recent progress in experimental and 
theoretical studies of Bose-Einstein condensates (BECs) 
[1] has made solitary matter waves physically relevant ob- 
jects. One-dimensional (ID) dark [2] and bright [3] soli- 
tons have been observed in recent experiments. On the 
other hand, optical solitons have a time-honored history 
as fundamental nonlinear excitations in optical fibers and 
waveguides (see, e.g., the recent reviews [4,5]). 

A very relevant generalization of this class of physical 
systems and of the solitary waves they can support con- 
cerns the case of multiple coupled components. There has 
recently been a considerable volume of work relevant to 
the properties of coupled BECs ranging from the ground 
state solutions [6,7] to the small- amplitude excitations 
[8] of the order parameters. Furthermore, the forma- 
tion of various structures including domain walls [9,10], 
bound dark-dark [11], dark-bright [12], dark- ant idark, 
dark-gray, bright-antidark and bright-gray soliton com- 
plexes [13], as well as spatially periodic states [14] was 
also predicted. On the other hand, experimental results 
have been reported for mixtures of different spin states of 
87 Rb [15] and mixed condensates [16,17]. It is relevant to 
mention the efforts into the realization of two-component 
BECs from different atomic species, such as 41 K- 87 Rb 
[18] and 7 Li- 133 Cs [19]. 

Typically, the relevant model for two coupled BECs 
involves two nonlinear ly coupled Gross-Pit aevskii (GP) 
equations. However, in experiments with a radio- 
frequency (or an electric field) coupling two separate hy- 
perfine states [15,20], the relevant model involves also a 
linear coupling between the wavefunctions. The govern- 
ing normalized equations are then of the form: 



#2t = 



A + V + a n |^i| 2 + a 12 |^2| 2 
A + V + a 12 |^i| 2 + a22^2| 2 



^i+a^ 2 , (1) 
+ (2) 



where V = V(r) is the relevant potential, typically con- 
sisting of a magnetic trap and/or an optical lattice [1,21], 



while -0j's represent the condensate wavefunctions. The 
intra- and inter-species interactions are characterized by 
the coefficients djj (j = 1,2) and au respectively, while 
a denotes the strength of the radio-frequency (or elec- 
tric field) coupling. Note that Eqs. (l)-(2), combining 
both linear and nonlinear couplings, occur in fiber optics 
as well: In that case, ^-'s are two coupled electric field 
envelopes of the same wavelength but of different polar- 
izations and the linear coupling is generated either by a 
twist applied to the fiber in the case of two linear polar- 
izations, or by an elliptic deformation of the fiber's core 
in the case of circular polarizations [22,23] (linear cou- 
pling is impossible when considering waves of different 
wavelengths). Another optical model, with only linear 
coupling between two modes, applies to nonlinear fiber 
couplers [24] or dual-core nonlinear fibers (see, e.g., [25]). 
In the context of BECs, this coupling has been recently 
examined for extended wave solutions in [26]. 

In the present work, we aim to study a unitary trans- 
formation in the context of Eqs. (l)-(2) that completely 
absorbs the linear coupling between the components into 
an oscillatory temporal dependence. We illustrate the 
value of this transformation in a two-fold way: on the 
one hand, we use it to understand the role of the linear 
coupling between the components as a means of creat- 
ing Rabi oscillations between the matter present in the 
two components (e.g., hyperfine states) [27] and their 
analog in power oscillations between polarizations in op- 
tical systems. On the other hand, we use it to con- 
struct exact time-periodic solutions of such linearly cou- 
pled nonlinear Schrodinger equations. This way, we iden- 
tify time-periodic Thomas-Fermi clouds and extended 
waves in one spatial dimension, as well as sloshing do- 
main walls and vortices in two dimensions. Since the 
transformation is exact only in the so-called Manakov 
case of an = &22 = a\2 [28], we investigate numeri- 
cally what happens to the constructed time-periodic so- 
lutions in cases (relevant to 87 Rb experiments) where 
an 7^ &22 7^ a\2. Finally, we generalize the transfer- 



mation to show that analogous constructions are feasible 
with a higher number of components. Our focus here 
is in illustrating the generality and applicability of the 
transformation to a very broad host of settings (multi- 
ple component, as well as higher dimensional cases, also 
in the presence of external potentials). We should note 
that special cases of this transformation have previously 
been identified in optics (in the integrable case of the two- 
component, one-dimensional Manakov solitons in the ab- 
sence of a potential [23]), as well as in a different format 
in BECs [in two- component, ID condensates in the pres- 
ence of a periodic potential [29], where solutions were 
constructed as unitary transformations of stationary so- 
lutions of the equations (6) below. There the nonsta- 
tionarity was introduced by the mixing of two stationary 
solutions that were not in phase.] 

Analytical Results. Assuming that an = <222 = g and 
ai2 = h, Eqs. (1,2) take the following matrix form: 

i^t ~ <*Pil> = + (^G^)ip + V(x)^, (3) 

where 

*-(?!;). *-(;;)■ «> 

Furthermore, in the special case h = g (hence G = gl, 
with / being the identity matrix), we may consider Eq. 
(3) as a homogeneous equation (LHS) with an inhomo- 
geneous part (RHS); in this case, the solution of the ho- 
mogeneous equation is given by 

/ tt/_l\ i -iaPt i ( cos(at) —ism(at) \ , 
^ = = e -»Vo = [ _ . s ^ t) cos( ^ } 'j*, 

(5) 

Using a variations of parameters approach, we substitute 
this in the equation for ip. This results in: 

#ot = -^A^o + (^oG^o)^o + ^(x)^ , (6) 

i.e., the same equation as (3), but without the electric 
field coupling terms. As mentioned above, large classes 
of exact stationary and nonstationary solutions to this 
equation were constructed in [9-14,29]. Those solutions 
can now be used to construct exact solutions to (3). Im- 
portantly, stable such solutions, are also relevant to the 
more general model Eqs. (l)-(2), and in particular to the 
BEC experiments with 87 Rb [16] or 23 Na [17]: in that 
case, the deviations of the values of the nonlinear coeffi- 
cients ajk from 1 are typically on the order of 3%. Such a 
small difference acts as a small perturbation and does not 
alter the stability of the solutions (see also the numerical 
results below). Note that, typically, these solutions will 
be nonstationary and, in particular, time-periodic. 

The success of this simple variation of parameters 
method can be phrased differently in physical terms: the 



equation (6) is invariant under unitary transformations, 
as it is an additive combination of the linear Schrodinger 
equation and a nonlinear term with a matrix G that com- 
mutes with any matrix. Then the extra linear coupling 
terms can be removed by using a time-dependent unitary 
transformation, which does not affect the other terms of 
the equation. Thus, the solutions of (3) can be thought 
of as rotating in "spin space" . 

We should note that more general nonlinear coupling 
matrices G commute with U. In particular: 

6-(H), ,7, 

also commutes with U. The relevant dynamical equa- 
tions, however, containing the term (^jG^o^o result in 
nonlinearities of the form \^i\ 2 ^2 and ^1^2 m the dy- 
namical equation for xpi and hence seem less physically 
relevant. 

The density of the different components rii = |^i| 2 
(i = 1,2) is given by 

m = cos 2 (at) I -001 1 2 + sin 2 (at) | ^02 1 2 

+ |^oi||^O2|sin(0 2 -0i)sin(2crt), (8) 

n 2 = cos 2 (at) I ^02 1 2 + sin 2 (at) | ^01 1 2 

+ l^oi 1 1^02| sin(0i -0 2 )sin(2af), (9) 

where ^01 = |^oi|e 26>1 , ^02 = |V^02 |e z6>2 . Note that if ^01 
and ^02 are time independent, then ri\ + ri2 does not de- 
pend on time. Thus it is a local constant of the motion 
(as opposed to the global constant J (n\ + ri2)dx which 
is conserved for any G). This is equivalent to the state- 
ment that the nonlinear term in (3) is invariant under the 
unitary transformation given by U(t). Notice also that 
the existence of stationary solutions of Eq. (3) requires 
the expressions (8)- (9) to be independent of t, which only 
happens when both ^10 and ^20 are stationary and equal. 
Furthermore, clearly the above unitary transformation il- 
lustrates that the linear coupling does not affect the inte- 
grable nature of the ID Manakov model (for V(x) = 0). 

Numerical Results. We now turn to the practical 
usefulness of the transformation i.e., constructing time- 
periodic solutions of Eq. (3) from stationary solutions 
of Eq. (6), as well as quasi-periodic solutions of the 
former from periodic ones of the latter (i.e., the trans- 
formation always inserts an additional frequency in the 
time-dependence of the solution). While one can follow 
this path also in the absence of the potential for the 
known, exact solutions of the Manakov model, we will 
focus herein on the case with the potential, which is more 
relevant to BECs [15-17,20]. 

Exact solutions in the presence of a potential are not 
often available in explicit form for Eq. (6). However, 
in the presence of the physically relevant, optical lat- 
tice potential [30] of the form V = — Vo sm 2 (mx) (in one 
spatial dimension), large classes of such stationary solu- 
tions can exist such as [14,31] ^01 = {y/Bi cos(mx) — 
i\JB\ + A\ sm(mx))e~ tUJlt , ^02 = (\^B2~cos(mx) — 
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i\JB2 + A2 sm(mx))e~ ZUJ2t , where (^1,2,^1,2,^1,2 are ap- 
propriate constants. These solutions then become exact, 
genuinely time-periodic solutions of Eq. (3) according 
to Eq. (5). Such a solution is given in the top panel 
of Fig. 1. The bottom panel of the figure illustrates 
a quasi-periodic solution of Eq. (3), constructed from 
a periodic one of the form ^01 = \fA\ sn(mx, k)e~ lUJlt , 
^02 = V~ A2 cn(mx, k)e~ tUJ2t in the elliptic function po- 
tential V = — Vo sn 2 (mx,fc), which degenerates into the 
OL one for k — > 0. 

However, typically such explicit solutions are not 
known, e.g., in the presence of a magnetic trap potential. 
Then one can use relaxational methods such as a Newton 
iteration or imaginary time integration to obtain station- 
ary solutions of Eq. (6) and exploit the time-dependence 
inherent in Eq. (5) to excite exact matter wave oscil- 
lations between the two components, producing exact, 
non-stationary solutions of Eq. (3). A one dimensional 
example of this strategy is shown in Fig. 2 for the ground 
state BEC in a magnetic trap potential [1,21] of the form 
V(x) = Vt 2 x 2 /2. ^2 is initialized at an exact stationary 
solution (in the presence of V(x), obtained via a Newton 
method), while ipi(x,0) = 0. However, an interesting 
question then concerns the potential persistence of the 
Rabi oscillations when the unitary transformation is no 
longer exact (i.e., for h ^ g). This is examined in the 
bottom set of panels in Fig. 2. Clearly, small deviations 
from the h = g limit (relevant to the BEC context where 
an ~ &22 « ai2) lead to persistence of the matter wave 
oscillations but with a quasi-periodic beating character. 
Beyond a critical threshold, however, the oscillations dis- 
appear and give rise to chaotic behavior. 

Naturally, the same idea for constructing nonstation- 
ary solutions can be carried over to two spatial dimen- 
sions. We illustrate the principle in the case of domain- 
walls (DWs) and vortices in two spatial dimensions (rele- 
vant results, but for a rotating trap were reported in [32]). 
As concerns the DW solutions, instead of using the more 
"standard" circular DWs (between a less repulsive com- 
ponent in the middle of the trap and a more repulsive 
one forming an outer shell) [6,7], we will use the recently 
proposed dipolar (i.e., rectilinear) and quadrupolar (i.e., 
cross-like) DWs of [33]. Once such a stationary state is 
reached for Eq. (6) (via imaginary time integration), us- 
ing it as initial condition in Eqs. (l)-(2) produces a time 
periodic solution of the latter. Such examples are shown 
in Fig. 3 for the different types of DWs in the case of a 
magnetic trap potential V(r) = (l/2)^ 2 r 2 (r 2 = x 2j ry 2 ). 
We should note that the DW can, in principle, exist if the 
immiscibility condition A = an (222 — a 2 2 < is satisfied. 
In particular, the quadrupolar DWs can persist even for 
A = (or g = h = 1) and hence the resulting "rotating 
propeller" solutions are exact time periodic solutions of 
Eq. (3). In the bottom panel of Fig. 3, a DW cross 
is shown for the practical case of an = 1.03, ai2 = 1 
and a22 = 0.97 (these values pertain to two different spin 
states of 87 Rb [16]). On the other hand, the dipolar DWs 
only persist if A < —0.061, which cannot be satisfied for 



the 87 Rb parameters. In the top panel of Fig. 3 such a 
dipolar DW is shown, in the limiting case A = —0.061 
(for g = 1 and h = 1.03); apparently, the time-periodic, 
dipolar DWs are only approximate solutions. Similarly 
to the ID results, if h ^ g (or, generally, A / 0), we 
have found that there is a maximum (minimum) criti- 
cal h c (A c ) up to which the Rabi oscillations persist [for 
the bulk of each component]: h c = 1.67 (A c = —1.78) 
for dipolar DWs and h c = 1.28 (A c = -0.64) for the 
quadrupolar DWs. 

Finally, one can use a similar construction for a pair of 
vortex structures with two components (see, e.g., Fig. 4). 
In this case, we initialize the imaginary time integration, 
in the absence of the linear coupling, with one compo- 
nent having a vortex centered at (5,0), while the other 
has a vortex at (—5,0). After the configuration relaxes 
to the stationary vortex pair solution of Eqs. (l)-(2), we 
again turn on the linear coupling and obtain a spiral ro- 
tation between the vortices resembling a spiral wave. In 
this case also, there is a critical h c = 1.32 (A c = —0.74) 
beyond which the regularity of Rabi oscillations is de- 
stroyed. In this case, the breakup leads to the formation 
of spiral patterns in the condensate. 

Conclusions and Generalizations. In this Letter, we 
have illustrated the possibility of coupled Bose-Einstein 
condensates to sustain exact periodic and quasi-periodic 
patterns in the presence of an experimentally realizable, 
linear coupling between the components. Similar results 
should be immediately applicable to linearly coupled op- 
tical systems close to the Manakov limit. A unitary 
transformation, commuting with the nonlinear kernel was 
identified as the source of such solutions and as a way of 
"factoring out" the linear coupling by means of time- 
dependent oscillatory behavior with a frequency equal to 
the strength of the linear coupling. We demonstrated 
the relevance of this transformation in constructing var- 
ious solutions in the presence of external potentials such 
as the optical lattice and the magnetic trap potential. 
We also illustrated the robustness of the mechanism in 
demonstrating that the phenomenon persists even for a 
wide range of parameter values (rather than only for the 
special yet experimentally relevant case of equal inter- 
and intra-species interaction for which it is exact). 

Finally, we would like to indicate that this mecha- 
nism is not restricted to the particular case of two lin- 
early coupled components, but in fact generalizes to 
higher numbers of components (e.g., 3 linearly coupled 
hyperfine states can also be realized in the context of 
BECs [26]). In this case, the term aPip in Eq. (3) 
should be substituted with where a is a symmet- 
ric n x n matrix (with zeros along the diagonal). The 
unitary transformation then becomes U(t) = e~ tta ; e.g., 
in the special case where n = 3 and all off-diagonal ele- 
ments are identical, then U(t) = [/%], with diagonal ele- 
ments fa = (l/3)(2e iat + e~ 2iat ) and off-diagonal ones 
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FIG. 1. Top panel: two spatial periods of the time-periodic 
oscillations of the density 712 for the trigonometric solu- 
tion in the OL potential. Here g = 1, at = 1, m = 1, 
lo\ — LO2 — \ra 2 + B x + B 2 , and V = 1, A x = 2, 
A2 = Vq — A\ — 1, B\ = 1, B2 — 2. Bottom panel: two spatial 
periods of a quasiperiodic-in-time oscillation of the density ri2 
for the elliptic function solution given in the text. Here m = 1, 
g = 1, a = 7r, (jj± = 7>m 2 (l + k 2 ) — A 2: 0J2 — |m 2 — A 2 . Also 
Vo = 1, A! = 2, A 2 = V + m 2 k 2 - Ai = 1. 
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FIG. 2. (Color Online) Top two panels: the two subplots 
show exact Rabi oscillations of a ground state BEC in the 
presence of a magnetic trap. The oscillation frequency is given 
by a = 0.8. The magnetic potential has a frequency Q = 0.1. 
The two subplots show the contour of the space-time evolu- 
tion of the density of each component. Middle two panels: 
evolution of the density at the center of the magnetic trap 
i.e., at x = for the two components. Bottom three pan- 
els: evolution at the center of the magnetic trap, for h ^ g, 
namely for h = 1.2 (top subplot), h = 1.5 (middle subplot) 
and h — 2 (bottom subplot). 




FIG. 3. Top panels: Contour plots of the density |^i| 2 for 
a dipolar DW (the density |V>2| 2 of the other species is com- 
plementary to |^i | 2 ) for t = (a), T/4 (b), T/2 (c), and 3T/4 
(d), with T n/a « 15.7 (a 0.2); Q 0.045, A -0.061. 
The pattern persists for long times with the species "inter- 
changing places". Bottom panels: Same as the top but for 
the quadrupolar DW with A = —9 x 10 -4 (the respective 
nonlinearity coefficients pertain to 87 Rb). The Rabi oscilla- 
tion of this sloshing DW gives the impression of a rotating 
propeller. 
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FIG. 4. Top panels: Contour plots of \ipi\ 2 for two cou- 
pled vortices, initially placed at x — ±5, for t = (a), T/4 
(b), T/2 (c), and 3T/4 (d), with T — 7r/a w 15.7 (a = 0.2); 
Q = 0.045, A = -9 x 10~ 4 ( 87 Rb). The vortices "interchange 
locations" (in a structure resembling a spiral wave). Bottom 
panels: Same as the top but for A = —3 (g = 1, h = 2). The 
configuration breaks up forming spiral patterns. 
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